Finite size effect of harmonic measure estimation in a DLA model: variable size of 

probe particles 
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A finite size effect in the probing of the harmonic measure in simulation of diffusion-limited 
aggregation (DLA) growth is investigated. We introduce variable size of the probe particles to 
estimate harmonic measure and extract fractal dimension of DLA clusters taking two limits, of 
vanishingly small probe particle size and of infinitely large size of a DLA cluster. We generate 1000 
DLA clusters consisting of 50 million particles each using off-lattice killing-free algorithm developed 
in the early work. The introduced method leads to an unprecedented accuracy in the estimation 
of the fractal dimension. We discuss the variation of the probability distribution function with the 
size of probing particles. 



I. INTRODUCTION 



There is a wealth of processes and systems ranging 
from the crystal growth and dendrites formation to bac- 
teria colonies growth and dielectric discharge patterns 
(see [1] for a review) that evolve in time in a similar 
manner following a generic process called the two dimen- 
sional aggregate growth. The specific model which is a 
most common representation of this process is the dif- 
fusion limited aggregation (DLA) Q and its generaliza- 
tion, dielectric breakdown model (DBM) ^ capturing 
most of the dynamic properties of the random aggre- 
gates evolution [ij . While analytical and numerical stud- 
ies have had impressively advanced our understanding of 
the DLA [iHl, several critical issues remain unresolved. 
One of the major controversies is related to the multi- 
scale and fractal nature of DLA. It has been a common 
belief that DLA, DBM and Laplacian growth [5j belong 
to the same universality class. Recently, however, this 
common wisdom was questioned in Ref. 0] , where it 
was argued that DLA and Laplacian growths have dif- 
ferent fractal dimensions. To resolve this intriguing issue 
a breakthrough in the precision of the data on fractal 
dimensions of both models is required. In our work we 
will address the algorithms and techniques that offer a 
dramatic improvement in determining the DLA fractal 
dimension D characterizing DLA clusters. 

The DLA/DBM models and growth algorithms are 
based on the aggregation of the particles randomly dif- 
fusing in the plane. A noticeable progress was made ow- 
ing to the technique introduced recently by Hastings and 
Levitov [8*1 where the electric field equations appearing in 
DBM were solved by iterative conformal mapping, and 
freely diffusing particles could reach any given site of the 
surface. The approach of ^ was further extended by 
Hastings @ where the models of the Laplacian Random 
Walk (LRW) similar to DLA, with the exception that 
the growth occurs only at the tip, were proposed. In the 



continuum limit this model transforms into the stochas- 
tic Loewner evolution (SLE) [l^, Yet, within the 
conformal mapping techniques fairly large clusters are to 
be built to insure a reasonable precision in determining 
the fractal dimensionality because the relative fluctua- 
tions decrease very slow with the cluster size. Namely, 
squared relative fluctuations 
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scale as iV-O-33 with the cluster size iV [H, [ll|. This is 
caused by large deviations in the cluster structure across 
the ensemble: each cluster has its own particular shape 
and preferred directions; these shapes vary in time and 
do not converge to any common "typical" structure. 

Random aggregate growth is fully characterized by the 
harmonic measure, i.e. by the probability for the surface 
to advance at some given segment. A common recipe for 
measuring this probability in computer experiments was 
to use the probe particles of the same size as particles 
that comprise the aggregate. To count the controversies 
concerning the issue of the universality of random aggre- 
gate growth, we undertake more precise analysis of DLA 
model, based on the probing particles of variable size. 
As it was shown this technique can signiflcantly in- 
crease the accuracy of the fractal dimension estimation. 
We build the ensemble of 1000 DLA clusters, 50 millions 
particles each, generated by the off-lattice killing-free al- 
gorithm [isj which we now modify to insure the free zone 
tracking. By using probe particles of different size S and 
then taking the limit D{S,N),S ^ we find fractal di- 
mension measured on the whole surface including also the 
regions that were inaccessible to the probe particles of the 
fixed size. Another important feature of our approach, is 
that the small particles are more sensitive to the fjords, 
while the large particles feel better the tips of the clus- 
ter. Thus, using the variable size-particles we measure 
not only the influence of the active growing zones (the 
hot part of the cluster) but also the effect of the "ac- 
complished" (frozen, or cold) parts of the cluster. As 
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a result, the large fluctuations of the fractal dimension 
D{S,N) are suppressed. Taking then the second limit 
D = D{N),N oo we find [H] that D converges to the 
value D — 1.7100(2) which is the order of magnitude 
improvement as compared to past simulations. 

The paper is organized as follows. In section [IT] we give 
details of the off-lattice killing-free algorithm we use to 
generate clusters and of the procedure of estimation of 
DLA fractal dimension. Section Hill contains brief discus- 
sion of the dependence of fractal dimension from the clus- 
ter size and section lTVl discusses walk of the center of clus- 
ter mass as cluster grows, and the two ways to calculate 
cluster radii - as distance from origin and as distance from 
center of mass. Section |V] describes the original method 
of fractal dimension estimation with the additional pa- 
rameter, the size of the probe particles. We discuss in 
section I VII how the probability distribution function for 
particle to stick cluster at the distance r varies with the 
size of the probe particles. Finally, in the section rvTIl we 
discuss our results and possible future developments. 

II. SIMULATIONS 
A. DLA algorithm 
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FIG. 1: Occupied lattice sites and sizes of free zones around 
each cell. The gray cell is the newly occupied one. Shaded 
cells are cells where free zone size will be recalculated. 



The original DLA algorithm |3| starts with placing a 
seed particle at the origin of a square lattice. Than at 
certain position, far away from the seed, a new particle 
is released, which wonders stochastically over the lattice 
until it runs into and sticks to the seed particle. If in 
the course of random walk a particle crosses the lattice 
boundary it gets removed (killed) . A new particle is gen- 
erated and process repeats. Successive aggregation of 
particles forms a DLA cluster. Several improvements to 
this original algorithm were developed in order to speed 
up the simulations [H, [l^ . Later it was found that large 
cluster exhibit anisotropy which reflects the symmetry of 
the lattice [17]. 

There exist also the off-lattice versions of the same al- 
gorithm. Namely, the particles are assumed to be balls 
that move freely in each direction. Making use the off- 
lattice algorithm one can easily speed up simulation s by 
utilizing improvements introduced in flal. [igI. [isl. [l9l [20|. 
in particular, the large step sizes, returns instead of 
killing, etc. 

Another algorithm for generating the DLA-like cluster 
is the well-known DBM model Q • Instead of tracing the 
motion of particles, one solves a Laplace equation and 
calculates a probability for a particle to hit the surface 
at the given site. Adding particles successively with the 
calculated probability, one generates a DBM cluster. 

B. Off- lattice killing- free algorithm 

We utilize the realization of the off-lattice killing-free 
algorithm described in fis'l, in which we upgraded the 



free-zone calculation. An essential feature of our algo- 
rithm, that it is a memory saving one; the memory or- 
ganization is similar to that used in the Ball and Brady 
algorithm [isj . The modifications we have introduced are 
as follows: (i) we employ the large walk steps; (ii) we use 
a killing-free rule for exact evaluating the probability for 
the particle to return at the birth circle; (iii) we use only 
two layers in the memory hierarchy; (iv) we use a recur- 
sive algorithm for the free zone tracking; and (v) the size 
of free zone is calculated precisely for the particles mov- 
ing near the cluster boundary. Our further improvement 
of the walk procedure [item (iv) above] is that instead of 
performing the search of free cells each time particle is 
moving, we pass all the cells and recalculate the free-zone 
sizes after the addition of each particle. 

The killing-free algorithm consists of the following 
steps: 1) Placing a seed particle at the origin (0,0). 2) 
Generating a new particle at a random position on i?^. 
3) Moving the new particle by fixed steps in random di- 
rections. 4) Aggregating a particle to the cluster in case 
it collides with either the seed or the cluster itself. 5) 
Generating a new particle position on Rf, with the prob- 
ability 



2tt — 2x cos S(p + 1' 

where x = r/Rh and 6(j) is the change in particle angle, 
in case where the particle falls into a dead circle Rd, 
Rd> Rb- 6) Repeating the cycle from the step 2. 

Since particles are moving randomly in the plane and 
assuming that the step size is negligible as compared to 
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the dimension of the free zones around a particle (un- 
less it gets close to the cluster), one could speed up the 
computational process by increasing the step size. This 
requires reliable and efficient identification of such "clus- 
ter free" regions. The common approach is to use the 
hierarchical memory model and store particles into cells, 
which will comprise the first layer. The second layer con- 
sists of cells which are several times larger than those at 
the first layer, and so on. Starting with the largest-scale 
cells layer, one checks their occupancy until the unoccu- 
pied cell is found. If the process is successful one gets 
the size of free zone at ones. Unfortunately, for the large 
enough clusters this approach is memory consuming and 
inefficient, because of the high degree of the granularity. 

We use a double layer memory organization: Particles 
are stored in the square cells, the latter also describe the 
size of free area around them. This information is up- 
dated during the attaching a new particle. Let i,j be the 
indexes of some cell on the lattice and Xn , Un be the in- 
dexes of some occupied cell. Then the free zone around 
the i,j cell is dij = min„TOaa;(|z — a;„|, \j — yn\) where 
the min is taken over all occupied cells. Figure [1] shows 
the lattice where occupied cells are black and numbers 
denote the free area sizes around each empty cell. As a 
new cell gets occupied (marked grey in the Figure) one 
has to check all the cells around it and recalculate dij. 
Shaded cells on the picture are those where the distances 
are going to be changed. One can see that the pertur- 
bation propagates continuously - all the shaded cells are 
connected with each other. This continuous propagation 
holds for every configuration of cells. 

In order to track down all the cells where the distances 
have changed, one should traverse first all the adjacent 
cells around the new one (grey cell in the Figure), then 
the next nearest cells that have changed the neighbors 
and so on. Since this process is unrestricted in time 
and can continue indefinitely (corresponding to the up- 
ward propagation in the figure), we stop the process after 
having spanned the distance of 15 and assume that cells 
where distance is not specified have free zone sizes of 15 
as well. 

For large clusters the cavities between the branches be- 
come bigger than 15 cell sizes. In order to overcome this 
limitation we use the two layer memory model each hav- 
ing free zone size counters updated by the rules we have 
described above. For clusters of size 50 million particles 
we choose first layer cell size of 32 and second layer - 128, 
and the particle radius is unity. 

On the close distance to the cluster we found useful 
to find precise distance to the cluster by iterating over 
all adjacent particles. This step is only required when 
current cell is marked as occupied or is adjacent to the 
occupied one and big step based on coarse grain informa- 
tion given by free/occupied cells is not enough. 

All these improvements of DLA realization algorithm 
enable us to generate each cluster with 50 million parti- 
cles in about 3 hours on 3Ghz Pentium 4 with 2 GB of 
RAM. 



C. Fractal dimension estimation 

The fractal dimension of a random aggregate is the 
quantity D defined by the scaling relation 

R(xN^/°, (2) 

where i? is a linear measure of the cluster consisting of 
N particles. Several possible choices of the characteristic 
length R are listed in Table |T1 



1. Ensemble averaging 

Let M denote the number of clusters we have generated 
and rM{k) be the position of the N-th particle in the A:-th 
aggregate. Then the deposition radius is 

1 

i^dcp(A^) = (rJv)-^5]r^.(fc). 

fe=i 

Other quantities we are interested in are the gyration 
radius, i?gyr, the root mean squared (RMS) radius i?2, 
and the penetration depth f , with the definitions given 
in the upper part of the Table U 

In order to understand the behavior of the fractal di- 
mension in the limit of infinite clusters one can check the 
dependence of D versus N. We extract D[Nq) as the 
result of fitting some R to the scaling relation ^ when 
iV varies from to A^o or in other words using clusters of 
size less than Nq. 



2. Harmonic measure averaging 

Harmonic measure is defined on the surface of the clus- 
ter as the probability of the DLA cluster to grow at the 
given point. One can define the deposition radius as 
^dcp ~ J '''d.q, where the integral is taken over the clus- 
ter surface and dq is the growth probability at point r. 
Using this harmonic measure we calculate fractal dimen- 
sion of fc-th cluster Dk{N), and perform averaging over 
the ensemble of clusters, D{N) = jyX^feli Dk{N). 

In simulations we estimate harmonic measure averages 
in the following way. We use probe particles which move 
due to the convention rules but do not stick to the clus- 
ter. Instead we track the position where i-th particle 
hits the cluster. Then deposition radius i?^™ estimated 

as average R^^ = Y^^iZi" ■ Calculating harmonic- 
measure averages is very time consuming procedure. The 
number of probe particles Np used is chosen dynamically 
so that relative error of i?^™ is less than 0.001. See the 
lower part of Table U for the definitions used. 
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Ensemble averages 
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TABLE L Estimates of tfie fractal dimension D extracted 
with the fit A'' oc R^ to the dependence of the various lengths 
with ensemble size of 1000 clusters each with 5 ■ 10^ particles. 
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FIG. 2: Fractal dimension as function of A'^ measured with 
i?dep and harmonic measure averaging over 100 clusters of 
size 5 • lO'^. Stars represent estimation in logarithmic scale, 
and boxes - in linear scale. In the inset; dependence of relative 
squared fluctuations of D with N . 



III. OSCILLATIONS OF FRACTAL DIMENSION 
WITH THE CLUSTER SIZE AND WEAK 
SELF- AVERAGING 

One can expect that the fractal dimension can be found 
as the limit D{N), as ^ 00. Unfortunately, this sim- 
ple approach fails because of the non-monotonic D(N) 
dependence with the rather large variations in D{N) [13|. 

We measure fractal dimension for each cluster using 
R'^^p and then we calculate the ensemble average of frac- 
tal dimension and its variation. There are two ways 
of extracting D from the R{N) dependence: either via 
choosing points in iJ^ep curve to be uniform on either 
the linear scale, or on the logarithmic one. It seems that 
both procedures result in the same average value of D, 
but exhibit different behaviors of the variances, see Fig- 



ure [S] If points are distributed uniformly in the logarith- 
mic scale, then jFp © changes as iV-0.33 found 
previously in our paper [3]. In contrast, D measured by 
the points uniformly distributed on the linear scale shows 
larger fluctuations without any visible decay with N. 

One can try to understand this behavior as follows. 
Varying points distribution we, accordingly, change their 
respective weights. If points are uniform on the linear 
scale, then D is more sensitive to the behavior of the 
bigger clusters, and the error when determining D does 
not change with the size of the cluster very much. One 
can say that D exhibits the lack of self-averaging in this 
case. On the contrary, the uniformity on the log scale 
gives rise to a weak self- averaging of D. 

The two procedures both yield very close values of D 
as one can check in Figure [H therefore for the rest of this 
paper we will use the points uniformity on the linear scale 
because this choice requires less CPU time for calculating 
the harmonic measure. 



IV. CENTER OF MASS FLUCTUATIONS 

Each realization of a DLA cluster is a unique object: its 
branch structure varies very much from sample to sam- 
ple. This causes large fluctuations in cluster properties 
and reflects in the behavior of the center-of-mass position 
of a cluster. At any moment of time (which is propor- 
tional to the number of particles in the cluster TV) there 
are some preferred directions of growth. These directions 
vary in time as illustrated in Figure [31 where we present 
the variations of the center-of-mass positions in plane co- 
ordinates {X, Y) for five different clusters. The center-of- 
mass positions (CoM) for each cluster are marked with 
their own symbol (circles, rhombi, triangles, stars, and 
boxes). Each mark placed after 2 • 10^ particles were 
added, so the history of the center-of-mass variation rep- 
resents the time from 10^ to 5 • 10^. One sees from Fig- 
ure[21 that for some clusters (marked with stars), position 
of the CoM is mainly rotated around the origin, while for 
some of them (marked with boxes) CoM is diffused from 
the origin in the given time interval. 

On average, clusters are quite uniform when looking 
for the average displacement of CoM. We estimate the 
average angle (0) = -0.01504998(±0.05535) of CoM (in 
the interval [— 7r;7r]) to be close to zero, while the angle 
variance is approximately a/ {(j)"^) — (0)^ = 1.76. This 
is very close to the value one can find assuming that 
is uniformly distributed in given interval: (0) = and 
v/(0^)-(<^)2 = VV3«1.81. 

There have been several attempts in the literature to 
relate the CoM distance from the origin Rm to the ex- 
pression defined as 



Rhiin) = 
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FIG. 3: (Color online) Center of mass position of five different 
clusters growing from 10^ to 5- 10^ particles, with the interval 
of 2 ■ 10^ See discussion for details. 
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FIG. 4; (Color online) Dependence of center of mass position 
Rm{N) with the cluster size N, as calculated averaging over 
1000 cluster with 5 • 10^ particles. Dashed line is a fit to 
Rm{N) and dotted line have a slope 1/2. 



where K is the number of clusters. We, however, find 
no reason for such a connection, because {X, Y) changes 
rather randomly. One can think of it as of a random 
walk in the plane, reflecting the competition between the 
growing activity of the branches. One thus expects the 
corresponding exponent to be equal to d = 2 (and not 
to D 1.71) in the limit of the infinite cluster size. 
Yet, the rare events such as the active branch growth 
effectively lower this value. We have found, that at large 
cluster sizes the generated exponent may be estimated as 
1.8558(13) which is notably larger that D « 1.71. The 
dependence of Rm{N) is plotted in Figure 2] with solid 
black line, the dash line is the fit, and dots line corre- 
sponds to the slope equal to 1/2 (reverse exponent). 



r to seed r to center-of-mass 

fidep 1.71111(59) 1.71112(60) 

7?2 1.71155(29) 1.71149(54) 

Rg 1.71149(30) 1.71133(30) 



TABLE II: Fractal dimension estimated by the standard 
method using the stick distance r to the seed (second col- 
umn) and to the center-of-mass (third column) as estimated 
for the number of typical lengthes. 

One can calculate the radius of the cluster by either (i) 
choosing the origin (0,0) as a reference point, or (ii) by 
choosing a temporary CoM position in the {X, Y) plane 
as a reference point in the estimation of radii given in 
Table HI It turns out that the fractal dimension calcu- 
lated by means of either of the recipes is the same within 
the computational error as one can see comparing the 
second and third columns in the table |TT1 However, this 
effect shows the additional possible reason for high fluc- 
tuation rates in DLA cased by the unique random branch 
structure of each cluster. 



V. FRACTAL DIMENSION ESTIMATION: 
VARIABLE PROBE PARTICLES 

A common technique of enhancing the precision of 
computations is the noise reduction in DLA. In doing 
so, one selects only the most probable events by attach- 
ing the hit counters to each particle of the cluster. The 
new particle is added to the particular old particle when 
the old particle is hit with some prescribed number of 
times. Similarly, the same procedure can be used when 
calculating the harmonic measure averages. Tips are the 
most probable places to grow, and at the same time the 
fluctuations at the tips are more strongly associated with 
the drastic changes in the cluster geometry, which is re- 
flected in the competing growing of the branches, in the 
birth of the new branches, and in the mutual screening 
of branches. The internal (frozen) part of the cluster sur- 
face, which is screened in the fjords and carry the small 
part of the measure, may give some contribution to the 
harmonic averages due to its length which can be large. 
So, we suppose that the internal part of the DLA cluster 
may contribute to the average. To check on this point, 
we introduce the variable size of the probe particles. We 
choose the size of the particles which we used to build the 
cluster to be the unit of the length, and measure in this 
units the probe particles size, 6. Larger particles increase 
the growth probabilities at the tips, while smaller parti- 
cles will penetrate dipper in the fjords, and redistribute 
harmonic measure shifting the part of the growth weight 
from the tips to the fjords (see next Section for details). 

Figure [5] illustrates the variation of the harmonic mea- 
sure estimation with the probe particle size S. Here the 
intensity reflects the probability (in the log scale) for the 
particle to hit the given segment of the cluster. The probe 



FIG. 5: (Color online) Fragment of the cluster with A'^ — 10^ particles. Intensity of the segments is proportional to the logarithm 
of the hit probabilities, corresponding to the probe particles of size 5 = 0.1, 1, and 10, from left to right. Probabilities are larger 
near the tips and smaller inside the fjords. For lower S the bigger part of the cluster surface is touched by probe particles. 



particles of the size (5 = 10 hit only few particles at the 
tips and never go inside the fjords. Only this hottest part 
of the cluster contributes to the harmonic measure esti- 
mation. The particles of the size 6=1 touch the bigger 
part of the surface but still there are regions inaccessible 
to them. The smaller particles of the size 6 = 0.1 (recall 
that this size is less than that of the building particles 
of the size 1) touch visibly the larger part of the cluster 
surface, penetrating deeper into the fjords. 



Two features contribute to the fjord screening: the 
size of a probe particle and the bottleneck configuration. 
Even if a particle fits into the bottleneck, the probability 
for it to pass is very small if the size of the particle is 
comparable to the bottleneck width. In other words, the 
effective channel width is somewhat narrower than the 
bottleneck width. Using the variable size probe particles 
we suppress the influence of the bottleneck and can thus 
measure the lower probabilities at the fjord surface. At 
the first glance, our approach looks opposite to the noise 
reduction method, which cuts the parts of the surface 
with the low harmonic measure. In fact, this contrast is 
misleading, because noise reduction method affects the 
growth process of the cluster, and our method affects the 
measurement process of the harmonic measure, and leads 
to the higher precision of the estimation of harmonic av- 
erages. 

One can think of our approach as of the canonic mea- 
surement of the length of the coast: the smaller the ruler, 
the bigger the coast length is [21] . In our case, the smaller 
probe particles, the longer the reachable surface is, with 
the cut-off on the scale of the bulk particles. 



A. Number of particles on the surface 

By varying the size of probe particles we change the 
reachable effective cluster surface. The number of parti- 
cles on the surface reachable during the probing process 
is shown in Figure [H Solid line on the picture is the fit 
to the expression 



Nreach^Nsurf/{l + S/5or (3) 

with a = 0.91(1), So = 2.23(6) and Nsurf = 45697(187). 
Nsurf is a total number of surface particles in cluster, 
those reachable by probe particles with infinitesimal size, 
and clusters size fixed to 10^. Thus, in the measure- 
ments with 10^ probe particles one can touch only about 
Nsurf ~ 45700. The number of reachable particles for 
(5 = 1 is lower by the factor of 1.4 than for particles 
of size S = 0.1. The number of reachable particles is 
the complex function of the probe particle size and of 
the number of probe particles. The best way to fill the 
whole surface of the cluster with probes is to find the 
limit N probes ~^ CO- But this is very time consuming 
procedure. 



B. D in the limit of vanishing probe particle size 

The fractal dimension of the aggregate turns out to 
depend also on (5, and the effective fractal dimension is 
now a function of two variables, D{N,6). We are inter- 
ested in the double limit D = liniN^oo.s-^o D{N,6). Let 
us first find the limit of the vanishing S, which gives us 
the D{N) dependence: 
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FIG. 6: Number of particles on the surface reachable by probe 
particles of size S. Measured by 10® probe particles and clus- 
ters of size 10®. 



FIG. 8: Fractal dimension as function of A'^ for the limiting 
value of (5 = 0. It reach value D = 1.7100(2) (plotted with 
the dotted line) in the limit of large cluster size A*'. 
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for the set of 6 is shown in the inset of Figure [71 the top- 
most curve corresponds to ^ = 100, while the lowermost 
- to ^ = 0.1. 

Result of the data fit to Expr. ^ are shown in Fig- 
ure [H 

Surprisingly, in the limit of the vanishing size of probe 
particles fractal dimension become monotone function. 
As N goes to infinity D{N) reaches its asymptotic value 
with approximate behavior (fit using the last 6 data 
points) D{N) = £'-1-4.2 • 10^7^^ " which gives the high 
precision of fractal dimension estimation D = 1.7100(2). 



FIG. 7: Fractal dimension D{S, N) as a function of the probe 
particle size 5 for different cluster size. Inset: D{S,N) is 
shown as a function of A'^ for the set of 5 values. Lines are the 
guide for the eye. See text for details. 



i:»(A^) = limi:>(iV,(5) (4) 
with the help of the fit 

D{S,N) = D{N)+AS^. (5) 

We calculated harmonic measure averages with the set 
of probe particles sizes: S = 0.1, 0.3, 1, 3, 10, 30, and 
100. The resulting fractal dimension D{S, N) is shown in 
Figure[7]as function of 6 for N = 5, 10, 15, 50 millions 
of particles. Fractal dimension D{S, N) as a function of N 



C. Errors estimation 

Making use of the least square fitting log(i?J™) to 
logiV when calculating the fractal dimension of a sin- 
gle cluster Dk{N,d), one can estimate an error on its 
determining. As has been mentioned above, the clus- 
ters are unique, and this results in large fluctuations of 
the fractal dimension over ensemble. This means that 
the error in the D{N,6) averaged over the ensemble de- 
pends mostly on its variation, while the errors in each 
Dk{N,S) are several times smaller. The next stage of 
out calculation is constructing the nonlinear square fit- 
ting (NLSF) to the equation ([5]) where all the data points 
have a weight, iversely proportional to the error, assosi- 
ated with the point (instrumental error in terms of Origin 
program). Resulting D{N) curve is then fitted by the ex- 
pression D{N) = D + A/N^ and the error for D is again 
estimated with the help of NLSF. 
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VI. PROBABILITY DISTRIBUTION 
FUNCTION 

Harmonic-measure averages can be written in a form 
-Rdep = / Pir, N)rdr, where P{r, N) is the probabihty 
for a particle to stick to the cluster with TV particles at 
the distance r from the origin (or from center-of-mass 
as discussed in Section IIVI) . This probability acquires a 
Gaussian form [12, [H, |24| | upon averaging over the en- 
semble of clusters. It is not so if one calculates P{r, N) at 
the surface of a single cluster realization; rare events as- 
sociated with the particular shape contribute essentially 
to the probability function. It was proposed in fTs*! that 
fluctuations in P{r, N) are responsible for multiscaling 
issue. 

We measure P{r, N, S) for a fixed cluster size N — 
2 • 10'' and the size of the probe particles 5, and aver- 
age over the ens emble of 1000 clusters. The resulting 
function P{r, N) is shown in Figure [H On the scale of 
the figure, all lines practically coincide. Fluctuations of 
the measured probability around the Gaussian form are 
visible despite the large number of clusters we use for 
averaging. 

The effect of the variation in the probe particle size 
on the probability function may be investigated with the 
use of the difference, 

SP{r,N) ^ P{r,N,d) - P{r,N,6 ^0), (6) 

which we approximate with the difference taken for the 
smallest used probe particle size S = 0.1, SP{r,N) = 
P{r, N, 6)-P{r, N,S = 0.1). 

The variation 5P{r, N) demonstrates the sin-like 
shape, with zero variation for any S at the distance cor- 
responding to the average deposition radius -Rdcp- The 
finiteness of the probe particles leads to the increase of 
the probability for r > Rdcp and probability decreases 
at distances r < i?dcp (shown with arrow in Figures [9] 
and [TUl The amplitude of SP{r, N) varies with S as 
^0.78(1)^ So we can normalize functions SP{r, N; 5) with 
the multipher G(r, TV, 6) = 6P{r, TV, 5)^-° ™. The 

set of this functions is presented in Figure [TOl All of 
them completely coincide well. The only exception is 
for 5 — 0.3 but this is clearly due to approximation 
P{r,N,5 ^ 0) w P{r,N,5 = 0.1) used in expr ©. It 
seems, therefore, that finite size of the probe particles 
will case the systematic overestimation of the fractal di- 
mension. 



VII. DISCUSSION 

Utilizing the variable probe particles for the harmonic 
measure estimations allows analysis of fractal dimension 
as a function of two variables. Fixing cluster size and 
taking the limit of the vanishing size of a probe particle 
effectively smoothes out the effects of the sample fiuctua- 
tions as is seen from the monotonic variation of the fractal 




r 



FIG. 9: Probability to stick at distance r for different 5 aver- 
aged over 1000 clusters with size A'^ = 2 ■ 10^. 
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10000 20000 30000 40000 

r 



FIG. 10: (Color online). Normalized change of probability to 
stick at distance r for different 5 relative 5 — 0.1. Fixed point 
is at -Rdop. Clusters size TV = 2 ■ 10 



dimension with the growing cluster size. This contrasts 
favorably the conventional methods where the fractal di- 
mension can oscillate with the cluster size making the 
convergence of the results very slow. As a result using 
variable size particles, one can extract the value of the 
fractal dimension from D{S,N) with the unprecedented 
high precision: we have found that D ~ 1.7100(2) for the 
off-lattice DLA cluster in the plane. 

We like to stress the accuracy of the above estimations 
as compared to those obtained by the traditional meth- 
ods mentioned in Tables |T] and |TT1 Technically, all the 
estimations are compatible within the 2a accuracy and, 
as it seems, may be treated with the equal confidence. 
The major difference of the values estimated with the 
help of varying probe size, is that it is estimated not as 
a single random variable but as the asymptotic value of 
some regular function as shown in Figure [51 Our ap- 
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proach corresponds to the variation of the ruler scale in 
the conventional estimation of the fractal dimension [2]| . 

The major result of our work is that for the first time 
we have reached a quality of the data which can offer 
reliable quantitative answers. In particular, the obtained 
value of fractal dimension allows to rule out the d= 17/10 
hypothesis proposed by Hastings about 10 years ago j25| . 
Furthermore, our technique can be used for the analysis 
of fractal dimension as a function of the lattice symme- 
try (d=3,4,5,6,7,8 and off-lattice corresponds to infinity). 
The next interesting question to address is the distribu- 
tion of the harmonic measure near the tips and inside the 
fjords, where one can expect our approach to be most ef- 



fective. 
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